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FOREWORD 

This document is am annual report describing the results of effort by 
personnel of Lockheed Missiles & Space Compamy, Inc., Huntsville Research 
& Engineering Center, for the National Aeronautics amd Space Administration 
under Contract NASW-3281, "Manufacturing in Space: Fluid Dynamics Nu- 
merical Analysis," The contractual effort described in this document was 
performed during the year from September 1981 to September 1982, The 
NASA Technical Director for this contract is Dr, Robert F, Dressier, 
Manager, Advanced Technology Program, NASA Headquarters, Washington, 
DC. 
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INTRODUCTION 

The research program described in this annual report is a continua- 
tion of effort initiated in August 1979 to investigate natural convection in 
various materials processing experiment configurations under microgravity 
conditions simulating the orbiting space station environment. The investi- 
gation involves the use of the Lockheed developed General Interpolants 
Method (GIM) fluid dynamics computer code to numerically simulate the 
development of natural convection under various loads. 

Two basic tasks were performed during this year’s effort; 

1. Last year’s investigation of the effects of changing the magnitude 
and direction of the acceleration vector was extended to include 
additional acceleration change sequences. 

2, An investigation was made of the effect of enclosure shape on 
natural convection velocities. 

The results of these tasks were published as separate documents 
which are included as Parts 1 and 2 of this annual report. 
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ENCLOSURES FOR VARIOUS ORIENTATION’S 
OF THE ACCELERATION VECTOR 
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ABSTRACT 

Numerical computations were performed for natural convection in circu- 
lar enclosures under various conditions of acceleration. It was found that 
subcritical acceleration vectors applied in the direction of the temperature 
gradient will lead to an eventual state of rest regardless of the initial 
state of motion. Supercritical acceleration vectors will lead to the same 
steady state condition of motion regardless of the initial state of motion. 

I 

Convection velocities were computed for acceleration vectors at various 
angles to the initial temperature gradient. The results for Rayleign 
numbers of 1000 or less were found to closely follow Weinbaum's first order 
theory. Higher Rayleigh number results were shown to depart significantly 
from the first order theory. Supercritical behavior was confirmed for 
Rayleigh numbers greater than the known supercritical value of 9216. 

Response times were determined to provide an indication of the time required 
to change states of motion for the various cases considered. 
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To 

AT 

t 

t* 

V 

'^rr.ax 

'^maxQ 
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NOMENCLATURE 


Description 

cylinder diameter = 2 R 
gravity force 

horizontal g component 
vertical g component 
radial distance 
cylinder radius 


g3 AT d 

Rayleigh number = ^ 

horizontal Rayleigh number 
vertical Rayleigh number 
temperature 

initial mid-point temperature 

temperture difference across circular cylinder 
time 

response time 
velocity 

spatial maximum velocity 

Vmax zero horizontal gravity component 
''max " ''maxg 

rectangular coordinates (Fig. 1) 
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Symbol 


NOMENCLATURE (Concluded) 
Description 


a 

B 

0 

U 

V 

P 


thermal diffusivity 

volumetric coefficient of thermal expansion 

polar angle (Fig. 1) 

dynamic viscosity 

kinematic viscosity , y/p 

stream function 

density 

p at T = - 
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INTRODUCTION 

The Inveetigation reported herein is bsiically a continuation of a pre- 
vious investigation reported in Ref* 1. In the previous study, ve investi- 
gated natural convection in circular enclosures in which the acceleration 
vector changes magnitudes and directions according to various schemes. The 
interest in the effect of time-varying acceleration loads on natural convec- 
tion is due to current interest in space processing experiments in which 
small time-varying spacecraft accelerations in orbit will produce convective 
stirring in contained fluids. The previous study considered a circular 
enclosure with the acceleration load consisting of two components, a verti- 
cal component, g^, applied normal to the initial temperature gradient and 
a horizontal component, g^^, in either the positive or negative direction 
parallel to the initial temperature gradient. The separate components were 
applied in the following sequences: (1) the vertical component applied 

until steady state, followed by the horizontal component superimposed to 
produce a new steady state; (2) both components applied simultaneously until 
steady state; and (3) the horizontal component applied first followed by the 
vertical component. This previous investigarion showed that, for subcri- 
tical Rayleigh numbers, the imposition of a given set of vertical and hori- 
zontal components yields the same steady state convective flow field regard- 
less of the order in which the two components are applied. The current 
study extends the pre.eding study to include the following conditions: (1) 

the vert Ic 'll component applied until steady state, followed by removal of 
the vertical component and application of a horizontal compon«^ut, and (2) 
application of a constant magnitude acceleration vector at various angles to 
the vertical direction. 

As in the preceding study, the Lockheed-developed General Interpolants 
Method (GIM) code was used in the numerical computations. The computations 
were performed on the NASA-Langley CYBER 203 system. 
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PROBLEM FORMULATION AND NUMERICAL SIMULATION 


The problem invettigated ia that of two-dimena tonal natural convection 
within the circular cylinder encloaure ahown in Fig. 1. 



Fig. 1 - Geometry for Circular Cylinder Encloau'i’e 

The initial temperature diatribution vaa baaed on a linear horizontal gradi- 
ent in the poaitive x direction, with the boundary pointa held conatant in 
time: 


T(r,o) - T^ + (AT/2) x/R 

T(R,t) • T + (AT/2) cob 0 
o 

The gravitation loading conaiated of two componenta, a vertical component 
in the negative y direction and a horizontal component in either the 
poaitive or negative x direction. 
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The nuatsrical eimuletion wee baaed on a nodal point diatribution gener- 
ated by the GIM code geometry module by apecifying an array of 20 x 20 area 
clementa with 21 x 21 nodal points. The circle waa traated aa a four-aidfid 
figure, each aide being a quarter-circle arc. The circular area waa divided 
into generalized quadrilateral clementa with curvilinear aides, with the 
m^dal points located at the four corners of each element. The circular 
g some try is shown in Fig. 2 with the computational grid network super- 
imposed . 

For convenience in the nuueTical simulation, we assumed a cylinder 
radius R of 1 cm and a temperature difference AT of 100 C. The gravity com- 
ponents were changed accordingly to yield the proper Rayleigh number 
values. The fluid was assumed to have the thermophyuical properties of 
water and to behave as a Boussinesq fluid in its thermrl expansion charac- 
teristics. The thermophysical properties used in the numerical simulation 
are listed as follows: 


Property 

Value 

Viscosity, U 

1 centipoise 

Thermal conductivity, k 

0.00143 cal/cm-sec-C 

Density, p 

1 gm/cu cm 

Specific heat, Cp 

i cal/gm-C 

Thermal expansion coefficient, g 

2.07 X 10“^/C 
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Fig. ? - Geometry and Computatonal Grid for Numerical Simulation 
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RESULTS 

Computed results are sho%m in Figs. 3 and 4 for the cases in which a 
vertical acceleration is applied until steady state followed by removal of 
the vertical acceleration and application of horizontal accelerations of 
various magnitudes in both the positive and negative x directions (see Fig. 
1). Vertical Rayleigh numbers of 1000 and 5000 were considered, and their 
results are shown in Figs. 3 and 4, respectively. As indicated earlier in 
Ref. 2, a vertical Rayleigh number of 1000 is withiu the low Raylei^ number 
range for agreement with the first order theoretical results of Weinbaum 
(Ref. 3), while a vertical Rayleigh number of 5000 is somewhat out of that 
range. Horizontal accelerations corresponding to horizontal Rayleigh num- 
bers of +1000, +2000, +5000 and +10,000 were considered. Results for the 
posicive horizontal Rayleigh number of 10,000 are uot shown on Figs. 3 and 4 
because of its supercritical nature, as will be discussed lacer. Note in 
both Figs. 3 and 4 that the negative horizontal Rayleigh numbers result in a 
slowing down of the initial counterclockwise circulating flow, leading to a 
reversal in flow direction and followed by a continuous approach to a state 
of rest. The positive horizontal Rayleigh numbers, however, result in a 
continual decrease in flow velocity approaching zero without a flow 
reversal. A simplified analytical model described in the appendix basically 
confirms the validity of these trends. 

Decay vesponse times were calculated for the results in Figs. 3 and 4 
to provide a quantitative indication of the time required for the flow to 
decay down to a state of rest after removal of the vertical acceleration and 
imposition of the horizontal acceleration. For the positive Ra^ values, 
the decay time is defined as the lime required for the spatial maximum 
velocity to reach the fraction 1/e (.368) of the initial value. Because of 
the ^low re' H . sal that takes place for negative Ra^ values, the decay time 
for th'se cases is defined as the time required after passing through the 
minimum (maximum negative value) to reach the fraction 1/e of the minimum 
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Fig. 3 - Spatial Maximum Velocity Histories After a Vertical 
Rayleigh Acceleration of 1000 is Removed and Various 
Horizontal Accelerations are Applied. 
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Approaches zero 
exponential decay 
time constant, 1.43 
(dimensionless) 


-1000 

-2000 

-5000 

- 10,000 


Dimensionless Time, Vt/R 


Spatial Maximum Velocity Histories After a Vertical 
Rayleigh Acceleration of 5000 is Removed and Various 
Horizontal Accelerations Are Applied 
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value. The decay times for initial vertical Rayleigh numbers of 1000 and 
5000 are shown in Fig. 3. Note that the decay times for both Rayleigh 
numbers are essentially identical except in the range near a horizontal 
Rayleigh number of 3000. Decay times for negative Rayleigh numbers range 

O 

from approximately 0.23 to 0.43 R /v, and for positive Rayleigh numbers 
from approximately 0.1 to 1.6 R /v. According to Weinbaum's theory, in 
the absence of a vertical acceleration component and for horizontal Rayleigh 
numbers less than the critical value of 9216, any velocity perturbation will 
be damped out, and convective motion will not develop. This critical value 
also follows from the simplified analysis in the appendix. In horizontal 
Rayleigh numbers greater than the critical value, however, a small 
perturbation :^n velocity will grow, and convective motion will develop. 

Shown in Fig. 6 is a comparison of convective motion resulting from a 
supercritical positive horizontal Riyleigh number of 10,000 applied to 
various initial states of motion. The initial states resulted from vertical 
Rayleigh numbers of 0, 1000 and 3000 applied until steady state and then 
removed. (It was necessary to initialize the fluid flow with a small 
perturbaton for the zero vertical Rayleigh number case for convective motion 
to develop.) Note that regardless of the initial state of motion, the 
convective motion proceeds toward the same steady state condition. In each 
case, an overshoot in the spatial maximum velocity occurs, followed by a 
continual decrease to the steady state value. 

A comparison is made in Fig. 7 of convective motion resulting from the 
two supercritical horizontal Rayleigh numbers of 10,000 and 20,000. A 
vertical Rayleigh number of 3000 was used to initialize the state of motion 
in both cases. The 10,000 horizontal Rayleigh number case is the same as 
that shown in Fig. 6. The 10,000 Rayleigh number case shows an initial 
increase and subsequent continued decay to steady state. The 20,000 
Rayleigh number case displays a strongly oscillating, although damped, 
approach to a steady state spatial maximum velocity approximately twice the 
10,000 Rayleigh number case. These limited results suggest a proportional 
relationship between horizontal acceleration aid the resulting spatial 
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Dimensionless Decay Time, vt /R 


5 



Horizontal Rayleigh Number, 


Ra^ 


Fig. 5 - Decay Time for Flow Initialized by Vertical Acceleration, 

Followed by Removal of Vertical Acceleration" and Imposition 
of Horizontal Acceleration 
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Fig. 6 - Comparison of Spatial Maximum Velocity Histories After Various Initial 

Vertical Rayleigh Accelerations are Removed and Horizontal Supercritical 
Acceleration ® 10,000) Is Applied 
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in/tximiiiii convective velocities for supercritical cases. Such a relationship 
would certainly break down near the critical Rayleigh number, however, since 
the steady state velocity goes to aero for Rayleigh numbers less than the 
critical value. 

The results shown in Figs. 8-10 are for acceleration vectors of cons- 
tant magnitude applied at various angles to the vertical direction. The 
magnitudes of the acceleration vectors correspond to Rayleigh numbers of 
1000, 5000 and 10,000. 

The steady state results are compared in Fig. 11 to the low Rayleigh 
number (first order) theory of Weinbaum which considers only the vertical 
component. The Ra > 1000 results are very nearly identical to Weinbaum* s 
theory, thus indicating "low Rayleigh number" conditions. The Ra » 3000 
results show considerable skewness in the curve, indicating deviation from 
low Rayleigh conditions, but still goes to zero at 0 ■ 90 deg, indicating 
subcritical conditions. The Ra * 10,000 results are clearly out of the low 
Rayleigh number range, and indicate supercritical conditions, since non-zero 
steady state velocities occur at 4* ■ 90 deg. 

Transient response times were calculated for the results in Figs. 8, 9 
and 10 to provide an indication of the time required for development of con- 
vective flow. Since some of the results showed an overshoot in velocity 
followed by decay to steady state, the response time was defined as the time 
required to reach the fraction 1-1/e (.632) of the maximum velocity acheived 
either at steady state or during the overshoot prior to steady state, which- 
ever is greater. The resulting response times are shown in Fig. 12 as a 
function of the horizonta'. Rayleigh number component for the various abso- 
lute Rayleigh number values. The simplified analytical development outlined 
in the appendix indicates that response times should be horizontal Rayleigh 

number dependent. The results shown in Fig. 12 indicate a general variation 

2 

in response time from approximately 0.03 R /v for negative horizontal 

2 

Rayleigh numbers to approximately 0.20 R /v for positive horizontal 
Rayleigh numbers. Note that, except for some scatter in the 3000 Rayleigh 
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n 

Dimensionless Time, vt/R 

Fig. 8 - Spatial Maximum Velocity History Starting from Rest for Ra - 1000 
and Various Off-Vertlcal Angles 
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Fig. 10 - Spatial Maximum Velocity History Starting from Rest for 
Ra - 10,000 and Various Off-Vertical Angles 
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Fig. II - Variation of Steady State Spatial Maximum Velocity 
with Angle Varying Between Vertical and Horizontal 
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Horizontal Rayleigh Number, Ra 8ln4> 


Fig. 12 - Transient Response Times versus Horizontal Rayleigh Number 
for Various Absolute Rayleigh Numbers 


17 


LOCKHEED-HUNTSVILLE PESEARCH & ENGINEERING CENTER 


LMSC-HREC TR D867640 


number date, the reeulte over the range of Rayleigh numbere coneidered 
closely follow a general t.'erd indicated by the faired curve. 

Contour plots of temperature, absolute velocity and stream function are 
shown in Figs. 13 through 16 for the case of a 1 cm radius circular enclo- 
sure containing water with a temperature difference c.; 100 K across the 
horisontal diameter. The acceleration vector was varied to correspond to 
the noted Rayleigh number values. 

The temperature contours at steady state are presented in Fig. 13 for a 

vertical Rayleigh number of 1000 and horizontal Bayieigh numbers of C, ^1000 

and -fSOOO. The positive horizontal components are shown to noticeably 
increase the distortion in the isotherms as expected. The negative hori- 
zontal components decrease the distortion somewhat, but not as noticeably 
in the positive case. 

The velocity contours at steady state are shown in Fig. 14, again for a 

vertical Rayleigh number of 1000 and horizontal Rayleigh numbers of 0, ;|_l0u0 

and +3000. The exhibited flow appears to h(‘. nearly perfectly circular in 
all cases with the maximum velocity indicated roughly near the r ' 0.3 to 
0.6 radial position. The theoretical maximum occurs at r « R/ vr- 0.38 R. 
Some of the departure from circular symmetry in the velocity contours is 
probably due to the numerical noise in -computing such small velocities. The 
circularity of the flow is more clearly indicated by the streamline contours 
shown in Fig. 13. A nearly perfectly circular flow pattern ic shown in all 
cases. 

The development of convective flow from rest for the supercritical 
Rs^ - 10,000 case is illustrated by the contour plocs of temperature tnd 
streamlines in Fig. 16. Steady state was reached in this case without break 

up into multiple convective cells because of the low value of Ra„. 

H 
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Temp*ratar* Contours 



Fig. 13 - Temperature Contours at Steady State for Various Horizontal 
Rayleigh Numbers (Ra^ « 1000, Temperature Units • K) 
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CONCLUSIONS 

The following general conclusions may be drawn from this study: 

1. Subcritical accelerations applied in the direction of the initial 
temperature gradient will result in an eventual return to a state of 
rest, regardless of the intensity of the initial state of motion. 

2. Supercritical accelerations apolied in the direction of the initial 
temperature gradient will result in an approach to the same steady 
state convective motion regardless of the initial state of motion. 

3. The variation in steady state velocities with angle at which the 
acceleration vector is applied closely follows Weinbaum's first 
order theory for absolute Rayleigh numbers of 1000 or less. For 
Rayleigh numbers greater than 1000, the deviation from Weinbaum's 
theory becomes significant. 
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Appendix 

SIMPLIFIED ANALYTICAL MODEL FOR NATURAL CONVECTION 
IN A HORIZONTAL TWO-DIMENSIONAL CIRCULAR ENCLOSURE 

A simplified analytical model was developed to predict natural convection 
In two-dimensional circular enclosures. The purpose of the simplified model Is 
to provide a basis for verifying the GIM code numerical results, and to assist 
In defining time constants for the computed transient response. The model 
development is similar to that reported in Ref. A-1. The geometry Is Indicated 
in Fig. A-1. 



Fig. A-1 - Geometry for Circular Cylinder Enclosure 
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llie initial temperature distribution was based on a linear horizontal gradient 
in the positive x direction, with the boundary points held constant in time: 


T(r,o) - T + (AT/2) x/R 
o 

T(R,t) - T + (AT/2) COS0 
o 


( 1 ) 


We assume the usual circular flow pattern with the velocity profile given by 


V 



Vx 


( 2 ) 


with v ^^^ varying with time during the transient flow period. We take the 
counter clockwise direction as positive, since the convective flow will assume 
this direction for a positive vertical acceleration component. 

We also assuo^ the transient temperature distribution given by 

AT 

T « T + ^ (r/R) COS0 + T' (3) 

O L 

where T* is a perturbation term due to convective flow given by 


T' 


iJl 

2 


T’ [1 - (r/R)^] (r/R) sin0 
max 


(4) 


with T ^^ varying with time. Along the y axis (0 ■ 7r/2), the perturbed temper 
ature distribution follows the same form as the velocity distribution given by 
Eq.2. 
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AtV 


The fluid denslit/ distribution is based on a Bousslnesq variation with 
temperature: 


P - pQ [1 “ 6(T - T^)] 


( 5 ) 


The torque. on the fluid mass due to viscous shear stress on Che Interior 
cylinder surface Is given by 


d V 

IV-aitRyji 


r - R 


“ 67t Vs y V R 
^ max 


( 6 ) 


The torque due to gravity Mg Is found by: 


«G 



gjj y) ds 


(7) 


where the Integration Is carried out over the circular region. The vertical 
gravity component g^ Is always positive and acts In the negative y direction. 
The horizontal component, g^. may be either positive or negative according 
to its direction with respect to the x axis. Carrying out the integration of 
Eq.(7) yields 


2tt R 




p(g^ COS0 + gg sln0) r^ dr 


d0 


( 8 ) 


o o 


f (*v * 
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ORIGINAL PhTE 3 
OF POOR QUALITY 

The Instantaneous angular momentum of the fluid mass Is given by: 

P ■ r pv dS 

where, again, the Integration Is carried out over the circular region. Carry- 
ing out this integration yields 


P 


2tt R 

’ffW 

o o 


dr de 


-7~ ttR p V 
4 o max 


( 10 ) 


By balancing the rate of change of angular momentum to the net torque on 
the fluid mass, a differential equation is obtained for the variation of v 

uBa 

with time: 


dP/dt - Mg + My 


& 

4 



P V 
o max 


J p SAT R^ 
o o 



- 6it V3 UR 


V 

max 


( 11 ) 


or 


V 

max 


1 ov 
16 V3 R^ 


Ra„ + V3 


T' 


max 

AT 





V 

max 


( 12 ) 


This differential equation contains two dependent variables, v and T* . A 

WBX i&ax 

second equation is thus required for a solution. We used the unsteady conduction 
equation. 
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OF POOR QUAUr/ 

+ V • VT - aV^T (13) 

evaluated at the point of maxinum velocity and temperature perturbation along 
the y axis, r ■ 1/V~3 (x • 0, y - 1/ /J), for the required second equation. 
Finding the temperature gradient and Laplarlan at x - 0 and y <■ 1//T from 
Eqs.(3) and (4) and solving Eq.(13) yields: 


T’ - V - 12 “ 
max 2R max 2 

K 

A single differential equation in terms of v^^ 
Eq.(12) and combining with Eq.(14): 


T’ 


max 

is obtained by differentiating 


(14) 


5 + ^ (1 

max Pr max Fr 


9216^ ''max" p^2 


0 


(15) 


where the tilda denotes normalization of the velocity to the characteristic 

2 

velocity, v/R; time to the characteristic time, R /v; and temperature to the 
temperature difference, AT. The general solution of Eq.(15) is 


where 


m t m t 

V - C.e ^ + C-e ^ + i-r- 

max 1 2 384VT 


Ra. 


' 9216 ^ 


(16) 


, (2 Pr 

m, • -6 r 

1 Pr 




(17) 


.filllL+Jl 

Pr 


1 


Pr 




(2 Pr + 1) 


2 " 9216^ 


] 


(18) 
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<=i - 


, (2 Pr + 1) 
« P? ‘■ 


1 +^ - 8 


Pr 


(2 Pr +1)' 


(1 


■Si->lk 


Ra 


wix^O “ 38^77 


Prd - 


1 2. (T* ) + 24 "v ) 

16 V5 Pr “ “* ® ' ° 




(19) 


6 <1IL.±J1 

Pr 


..fT 


Pr 




(2 Pr + 1) 


(1 li-) 

2 9215 ^ 


(v__) - 

0 




- 384V3 


Pr(l - 


9216 


■*■ pj. Te“ ^^liax^O ■ ^^max^O 


12 




Pr 




(2 Pr + 1) 


_ (1 ii- ) 

2 9216 ' 


( 20 ) 


where the zero aubscrlpta denote Initial conditions. 

We are particularly concerned with the following sequence. First a vertically 
downward gravity force Is inposed until a steady state condition is reached. The 
vertical gravity force is then removed, and a horizontal gravity force is imposed 
until a new steady state condition Is reached. Considering the Instant of imposi- 
tion of the horizontal force as the starting time, the initial conditions are as 
follows: 
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( 21 ) 

( 22 ) 


where (Ra )^ Is the Rayleigh number corresponding to the Initially imposed 
V ” 

vertical gravity force. 


For this case. Eq.l6 becomes: 
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For Pr > 7 (water). Eq.(23) reduces to the following double exponential equa- 
tions for the noted Ra^^ values. 


e Ra ■ -5000 

n 


- X / rt ilo/ -2.76 t ^ -22.96 t. 

V /(v )« - (- 0.694 e + 1.69 e ) 

max ntax 0 


(24) 


e Rajj - +5000 


V 

max 


/(v ) 

max 


0 


(0.577 e 


-0.772 t 


+ 0.423 e 


-24.95 t 


) 


(25) 


These equations are shown plotted in Fig. A-2 compared with GIM code numerically 
computed results for initial vertical Rayleigh number. (Ra^,)Q. values of 1000 
and 5000. Note that the analytical curves predict a flow reversal prior to an 
eventual decay to a state of rest for the negative horizontal Rayleigh number 
cases (Ra^ “ -5000). The GIM code results show good agreement with the analytical 
results for Ra^ ■ -5000 for both the 1000 and 5000 (F^)q values. Good agree- 
ment is also shown for the ■ 5000. (Ra^) " +5000 case. * 1000. 

Ra^ - +5(X)0 cape does not exhibit the same degree of agreement as the other 
cases, but the exponential decay trend is roughly the r?fle for both the analytical 
and numerical results. 
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ABSTRACT 


A numerical analysis was performed to compare natural convection 
velocities in two-dimensional enclosures of various shape. The following 
shapes were investigated: circle, square, horizontal and upright 2x1 

aspect ratio rectangles, horizontal and upright half-circles, diamond 
(square oriented with diagonal vertical) and triangle (equilateral with 
horizontal base). In all cases, the length scale in the various dimen- 
sionless parameters, such as Rayleigh number, is defined as the diameter of 
the equal area circle. Natural convection velocities were calculated for 
Rayleigh numbers of 1000 and 5000 with the temperature difference taken to 
be across (a) the maximum horizontal dimension, (b) the median horizontal 
line (line through centroid) and (c) the horizontal distance such that the 
temperature gradient is the same for shapes of equal area. A Rayleigh 
number of 1000 is within the "low Rayleigh number" range for agreement with 
first order theory for circular enclosures. A Rayleigh number of 3000 is 
slightly out of this range. For the class of shapes including the square, 
upright half-circle and upright rectangle, the computed velocities were 
found to agree very closely with that of the equal area circle when the 
temperature difference is taken to be across the maximum horizontal 
dimension (condition (a)). The velocities for the horizontal rectangle and 
half circle were found to be approximately one-half that of the equal area 
circle for the same condition. Better overall agreement among all shapes 
was obtained by setting the temperature difference across a distance such 
that the temperature gradients were equal for shapes of equal area. 
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Ra 

AT 

t 
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V 
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NOMENCLATURE 

Description 

equal area circle diameter 
gravity force 
circle radius 

3 

Rayleigh number = 

va 

temperature difference 
time 

velocity 

thermal diffusivity 

volumetric coefficient of thermal expansion 

dynamic viscosity 

kinematic viscosity, U /P 

stream function 

density 


PRECEDING PAGE BLANK NOT niMED 


Vi 


LOCKHEED-HUNTSVILLE RESEARCH & ENGINEERING CENTER 


LMSC-HREC TR D867640 


1. INTRODUCTION 

A number of fluid mechanics experiments are to be performed aboard 
orbiting spacecraft, primarily to investigate the possible use of the 
near-zero gravity environment in various materials processing applications. 

A knowledge of the intensity of convective stirring due to residual accele- 
rations is required in order to properly plan and evaluate the experiments. 
Various shaped containers will be used in these experiments. Theoretical 
results have been obtained which yield exact predictions of natural convec- 
tion velocities in an idealized container, i.e., the two-dimensional cir- 
cular enclosure (Ref. 1). It would be helpful for estimation purposes if a 
reasonable means existed for the extrapolation of the circular enclosure 
results to more complex shapes. A previous study demonstrated extremely 
good agreement in computed natural convection velocities for circular and 
square enclosures of equal cross-sectional area (Ref. 2). Dressier (Ref. 3) 
made use of this noted agreement in circular and square enclosure results, 
in addition to the more general results documented herein, to analyze 
natural convection in the proposed Lai, Kroes and Wilcox crystal growth 
experiment to be performed infight on Spacelab 3. The purpose of this 
investigation is to develop reasonable extrapolation criteria by comparing 
numerically computed natural convection velocities for various two- 
dimensional enclosure shapes with the circular enclosure results. The 
Lockheed-developed General Interpolant Method (GIM) computer code (Ref. A) 
was used in the numerical computations. 
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2. METHOD OF APPROACH AND NUMERICAL SIMULATION 


The following seC of two-dimensional enclosure shapes was considereo: 
circle, square, 2x1 aspect ratio rectangle in both upright and horizontal 
orientations, half-circle in both upright and horizontal orientations, 
diamond (square oriented with diagonal vertical) and triangle (equilateral 
with horizontal base). The baseline condition was selected for equal area 
with the temperature difference set across the maximum horizontal dimen- 
sion. The initial temperature distribution was based on a uniform hori- 
zontal gradient with the boundary points held constant in time. The gravity 
vector wac considered to be constant in the downward direction. The fluid 
was assumed to behave as a boussinesq fluid in its thermal expansion 
characteristics. 

The numerical simulation was based on dividing the various enclosures 
into a computational grid consisting of a network of generalized quadri- 
lateral elements with curvilinear sides, with the nodal points located at 
the four corners of each element. Each of the enclosures considered in this 
study was treated as a generalized quadrilateral region divided by inter- 
polation into an array of 20 x 20 elements. For example, the circle was 
treated as a four-sided figure, each side being a quarter-circle arc. 

The triangle was treated by forming a parallelogram such that two 
adjacent sides and the diagonal formed the desired triangle. The entire 
parallelogram was divided into a computational grid, but the diagonal points 
were treated as boundary points. The parallelogram with the diagonal 
boundary thus formed two independent triangular regions. The computational 
grids for the enclosures considered in this study are sho%m in Figs. 1-5. 
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Ihe grids for Che square and diamond are Che same since Che diamond is 
simply a square with Che diagonal in Che uprighc posicion. The recCangle 
and half-«ircle grids are idencical in eicher Che uprighC or horisoncal 
posicions. 
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3. RESULTS 

The velocity histories for the various shaped enclosures are sdown in 
Figs. 6 and 7 for Rayleigh numbers of 1000 and SOOO. The various enclosures 
are all of equal area, with the dimensionless velocities and Raylei^ 
numbers defined based on a temperature difference ilT across the maximum 
horizontal dimension, and a length scale equal to the diameter d of the 
equal area circle. Note that, for a Rayleigh number of 1000 the steady 
state velocities of the equal area circle and square are very close to 
identical, with the upright rectangle and half-circle being within about 5% 
of the circle value. The horizontal rectangle and half-circle, and the 
triangle steady state velocities form another grouping of values approx- 
imately one-half the circle value. The diamond results appear approximately 
midway between the two extremes. Roughly the same trend of steady state 
values appear for a Rayleigh number of 3000. Overall, the steady state 
velocities for the various equal area enclosures are in agreement within a 
factor of approximately 2. 

Response times for the various enclosure results are shown in Fig. 8 
for both 1000 and 3000 Rayleigh number. Agreement within a factor of 
approximately 2 is shown for all enclosures. These response times are 
defined as the time required to reach the fraction 1 - 1/e (0.632) of the 
maximum or steady state velocities, whichever is greater, in Figs. 6 and 7. 

An attempt was made to find other bases for correlating the data to 
yield better agreement between the results for the various shaped enclo- 
sures. Including the base line correlation, described earlier, correlations 
were made for the following sets of conditions concerning the distance over 
which the temperature difference was taken: 

• Temperature difference across maximum horizontal dimension 
(base line correlation). 
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• Temperature difference acroaa median horisontal Line (line 
through centroid). Only the triangle and horizontal 
half-circle reaulta changed in thia correlation from the 
baaeline correlation. 

e Temperature difference taken acroaa a diatance auch that the 
temperature gradient ia the aamc for all ahapes. 

The ateady atate reaulta for the above aeta of correlationa are aum- 
marized in Pig. 9 for a Rayleigh number of 1000, and in Fig. 10 for a 
Rayleigh number of 3000. For each of the correlation aeta, the length scale 
in the Rayleigh number and dimenaionleaa velocity ia taken to be the diam- 
eter of the equal area circle. A comparison of the three correlation sets 
shows improved correlation for the triangle and horizontal halt-circle by 
taking the temperature difference across the median horizontal line rather 
than the maximum horizontal dimension. Generally better correlation is 
obtained by setting the temperature gradients equal. This is equivalent to 
taking the temperature difference across a distance equal to the diaoioter of 
the equal area circle. 

Computer generated streamline, absolute velocity end temperature con- 
tour plots for all of the enclosures are shown in Figs. 11 through 34 at 
steady state for both 1000 and 3000 Rayleigh numbers. 
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4. CONCLUSIONS 

Natural convection velocities within two-dimensional enclosures of 
various shape, for Rayleigh numbers up to at least 5000, may be estinw id 
with reasonable accuracy by considering the area to be equivalent to a 
circle of equal area. For the class of fi.gures including the square, 
upright half-circle, and upright 2x1 aspect ratio rectangle, excellent 
agreement is obtained by considering the temperature difference across the 
maximum horizontal dimension to be equal to that across the equal area 
circle horizontal diameter. The agreement of the two upright oblong shapes 
indicates probable agreement for any similar upright oblong shape of roughly 
the same aspect ratio. 

Die horizontal oblong shapes have natural convection velocities 
approximately one-half that estimated based on the equal area circle with 
the temperature difference across the maximum horizontal dimension taken to 
be equal to that across the equal area circle diameter. As with the upright 
shapes, this probably indicates a general pattern for similar oblong shapes 
of roughly the same aspect ratio. 

Better overall agreement among all shapes is obtained by setting the 
temperature difference across a distance equal to the diameter of the equal 
area circle, thus making the temperature gradients equal for any shape. 
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Fig. 1 - Computational Grid: Circle 



Fig. 2 - Computational Grid: Square and Diamond 
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Fig. 3 - Computational Grid: 2x1 Aspect 3.0 Rectangle 



Fig. 4 - Computational Grid: Half-Circle 
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Fig. 5 - Computational Grid: Triangle 
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Fig. 6 - Velocity Histories for Natural Convection in Two-Dimensional 
Enclosures of Various Shape and Raleigh Number of 1000 
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Fig. 9 - Comparison of Steady State Velocities for Natural Convection 
in Two-Dimensional Enclosures of Various Shape and Rayleigh 
Number of 1000 
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Fig. 10 - Comparison of Steady State Velocities for Natural Convection 
in Two-Dimensional Enclosures of Various Shape and Rayleigh 
Number of 5000 
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Fig. 13 - Temperature Coutour Plots for Circular Enclosure 
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Fig. 15 - Absolute Velocity Contour Plots for Square Enclosure 
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Fig. 16 - Temperature Contour Plots for Square Enclosure 
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Fig. 
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Fig. 19 
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Absolute Velocity Contour Plots for Horizontal 2x1 
Aspect Ratio Rectangular Enclosure 



Contour Plots for Horizontal 2x1 Aspect Ratio Rectangular Enclosure 




Absolute Velocity Contour Plots for Upright Half-KIircle Enclosure 
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Fig. 25 - Temperature Contour Plots for Upright Half-Circle Enclosure 
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Absolute Velocity Contour Plots for Horizontal Half-Circle Enclosure 
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Fig. 30 - Absolute Velocity Contour Plots for Diamond Enclosure 
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Fig. 31 - Teaperature Contour Plots for Dianond Enclosure 
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Fig. 34 - Temperature: Contour Plots for Triangular Enclosure 


